Propagating spinors on a tetrahedral spacetime lattice 
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We derive a discrete path integral for massless fermions on a hypercubic spacetime lattice with 
null faces. The amplitude for a path with N steps and B bends is ±(l/2) JV (i/v / 3) s . 



In this article, we discuss a four-dimensional variant of 
'Feynman's Checkerboard' problem |lj — a construction 
of the propagator for the Dirac equation in 1+1 dimen- 
sions via a sum-over-paths method, where the paths tra- 
verse a null lattice. In that problem, the amplitude for 
a given zig-zag path is just (iem) R , where e is the time 
duration of a lattice step, m is the particle mass, and R is 
the number of direction reversals (and we use units with 
H = c = 1). If the mass vanishes the particle moves only 
to the right or the left with the speed of light, and these 
two motions correspond to the two chiralities for a Dirac 
spinor in 1+1 dimensions. 

Feynman's checkerboard path integral is striking in 
its simplicity, and intriguing because it accounts for rel- 
ativistic propagation and Dirac matrices with nothing 
more than a simple factor of i associated with a geomet- 
ric property — a bend — in a piecewise lightlike path. It 
hints at the possibility that a simple, discrete dynamics 
could underly the usual continuum description of rela- 
tivistic quantum field theory. A longstanding question 
has been whether or not a path integral could be found 
in 3 + 1 dimensions preserving these surprising features. 
Here we show that one can indeed come very close. 

Most work on lattice formulations of spinor propaga- 
tion is directed at lattice field theory calculations, and 
thus involves Grassmann variables and "path" integrals 
over field configurations in a spacetime of Euclidean sig- 
nature 0]. Here instead we seek a bosonic (i.e. non- 
Grassmanian) formulation in terms of a sum over particle 
paths on a Minkowski signature lattice. Such formula- 
tions have been studied previously in 3, 4] and references 
therein. What is new here is the interesting structure of 
the lattice employed, and the fact that it allows for a 
particularly simple rule for the amplitudes. 

We focus on the massless case, which in 3+1 dimen- 
sions is far more interesting than in 1+1 dimensions. The 
decoupled right and left chiralities are then described by 
two-component spinors satisfying the Weyl equation, 

o-»d^ = 0, (1) 

where — (1,+cr) with a the Pauli matrices. (The + 
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sign corresponds to right handed spinors.) 

The lattice we use here is topologically hypercubical 
but "tilted on its corner" . A quantum network based on 
this lattice was called "hyperdiamond" and previously 
studied by Finkelstein and collaborators^ (see also 0)- 
We use this name here for the lattice itself. 

A diagonal of the hypercube defines our time axis. The 
edges from one vertex lie in the direction of the four 
spacetime vectors 



nf = (1, afii). 



(2) 



The four spatial unit vectors n, point to the vertices of 
a tetrahedron, and the step speed a is for the moment 
unspecified. The hi sum to zero, and the inner product 
or angle between any distinct two is the same. Hence for 
i j we have hi ■ hj = —1/3, and the angle is equal to 
cos-^-l/S) re 109°. 

The spatial lattice at one time is a face-centered cu- 
bic (fee) lattice. A way to see this is to begin with the 
tetrahedron of points that lies at one time step to the 
future of a given spacetime point p. The four dimen- 
sional lattice has translation symmetries that map any 
point to any other point, and the spatial lattice at one 
time must share this property. Hence it can be grown 
from this tetrahedral seed by translation along the edges 
of the tetrahedron, which produces the fee lattice shown 
in Fig. HJ 




FIG. 1: Face-centered cubic lattice of points at one time step. 
The tetrahedron (dotted lines) is comprised of the four points 
reached from the center of the small cube in one time step 
(dashed lines). The continuum sphere of light is enclosed by 
and tangent to the tetrahedron. The distance from the center 
to a tetrahedron vertex is three times the radius of the sphere. 
The step length a and cube edge length L are shown. 

Evolving the spatial lattice one time step to the future 
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amounts to shifting it along the displacement from the 
center of one tetrahedron to one of its vertices, yielding 
a distinct but equivalent spatial lattice. After four such 
steps the original spatial lattice is recovered. 

It might seem natural to choose the step speed a = 1 , 
so that the links of the hyperdiamond would be null as 
envisaged in 0. However, in this case the lattice prop- 
agator would fail to converge at all in the continuum 
limit. The reason is that such a spacetime lattice vio- 
lates the well-known "Courant condition" for stability: 
the discrete region of causal influence must contain the 
continuum one. 

To marginally satisfy the Courant condition, the poly- 
hedral cone formed by the four hypersurfaces spanned by 
three of the n^'s must be tangent to the continuum light 
cone. That is, the faces of the hypercube must be null. 
Thus each of these hyperplanes must contain one and 
only one null direction. By symmetry this null direction 
must coincide with the sum of the three link vectors, e.g. 
A" 1 = (3, a{hi + h 2 + h 3 )). The Minkowski norm of 
is 9 — a 2 , hence if it is to be null we must choose a = 3. 
Moreover, N^n^ = 3 - (a 2 /3), so if a = 3 the null 
vector is orthogonal to all vectors in the hyperplane, 
confirming that the hyperplane is indeed null. 

Now consider the tetrahedral quartet of unit vectors 
hi. The sums ^ hf and J2 i n"n£ (with 'n°' denoting the 
'a'component of hi) are invariant under the symmetries 
of the tetrahedron, hence the first sum must vanish and 
the second sum must be proportional to the Euclidean 
metric 5 ah . Since the trace is equal to four this yields the 
relation 

J2^ = ^ ab - (3) 

i 

Using these identities and the definition @ of the 4- 
vectors nf, the matrix 4- vector = (l,c?) can be ex- 
pressed as 

^4eK i+ oM<- (4) 

In the special case a = 3, for which the polyhedral light 
cone conists of null hyperplanes, this becomes just 

i 

where 

P t = \{l + hi -a) (6) 

is the projector for spin up in the direction hi. Using 
the identity J5J the Weyl equation for right-handed 
spinors takes the form 

^^$ = 0. (7) 



We consider the hyperdiamond lattice with step vec- 
tors erij scaled by the step size e. With the partial deriva- 
tives replaced by finite differences, 

en t " l d^{x) = ^{x)^^{x-tn l ), (8) 

the Wcyl equation J7J yields the one-step evolution pre- 
scription for determining ^(x) on the lattice from the 
values — ern) at the immediately preceding points, 

%) = ^P,%-e4 (9) 

i 

The finite difference equation is different from any 
typically used in lattice field theory calculations. For a 
plane wave solution of the form ^(x 1 ) = exp^ik^x^^o 
(with is a constant spinor) it implies the dispersion 
equation 

i 

Expanding in e and using equation (J5J , we find at first or- 
der in e the standard Weyl equation for momentum eigen- 
states, a^k^o = 0. The solutions obey the relativistic 
dispersion relation, fc M /c M = 0, and the corresponding zero 
frequency solution has vanishing wave vector. 

Fermion doubling does not occur. That is, there are 
no additional zero- frequency solutions to Eq. (|10[l. as we 
argue later. The Nielsen-Ninomiya theorem [2| is pre- 
sumably evaded since this finite difference equation does 
not satisfy the hermiticity condition. The attendant lack 
of unitarity would be a problem for a lattice field theory 
application, but it does not present a problem for our 
purpose, which is just to extract a representation of the 
propagator as the continuum limit of a sum over paths 
on the lattice. 

Iterating Eq. (JSJ backwards in time, we see that the 
retarded propagator between two points can be written 
as a sum over multi-step paths involving moves in the 
directions h^ . . . hi N at step speed 3c, with the amplitude 
for such a path given by the operator 

~^[PiN ■ ■ ■ Pii ■ (11) 

The propagator is then the sum of these operators over 
all paths that connect the points. In terms of the unit 
eigenspinor \i) of Pj = the amplitude l(TT|l becomes 

<«iv |«V-i> <«jv-i | • • • |*2><*2|*i><*a|- (12) 

Note that for fixed initial and final points the number of 
each of the four step types is fixed, so one sums only over 
the order in which the steps are taken. 

The amplitude i|12|) is independent of the choice of 
phases for the spinors, so we are free to adjust those 
phases to produce a particularly nice result. The unit 
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spinor corresponding to a unit vector with spherical an- 
gles (0,(j)) is (cos(0/2), sin(0/2) exp(«0)) times an arbi- 
trary overall phase. The spinors \i), which correspond to 
the four unit vectors pointing to the vertices of a tetra- 
hedron, can be taken as 



|2) 



Jlp2 



1 



\V2 

i 4 > = 7sU^ >'"" 



If we choose 0i = — tt/2 and 02 = "03 — "04 = 0, then 
the inner products between the various spinors become 
identical up to a sign that depends on the order: 

(1|2,3,4) = (2|3) = <3|4) = (4|2) = -^. (14) 

Thus the amplitude for an N step path with B bends is 

B 

(15) 



± 



1 / i Y 



2 N \V3j ' 



To eliminate the operator | zjv) {* 1 1 from the amplitude 
(|12fl , we considered here the matrix element of the prop- 
agator with respect to fixed inital and final spinors \io) 
and |zjv+i) selected from the set ff 1 31) . These specify di- 
rections of arrival at the initial point and departure from 
the final point, and the amplitude 1|15|1 includes any con- 
tributions from initial and final bends. Unfortunately 
we have not been able to find a simple way to specify 
the overall sign other than by reference to the positive 
two-step orders given in (1141) . 

For left handed Weyl spinors the amplitude for a step 
is given by the orthogonal spin projector relative to the 
right handed case. Hence we can use the charge conju- 
gates of the four spinors 1)13(1 . with the result that the 
imaginary unit i in the amplitude (|15f) is replaced by — i. 
The effect of a mass m can be included by allowing for 
chirality flips between right and left handed spinor prop- 
agation at each time step, with an associated amplitude 
factor iem Q , as on Feynman's checkerboard. 

We now evaluate the sum over paths for the lattice 
propagator K e (Ax) for a spacetime displacement Ace and 
demonstrate that it reproduces the continuum propaga- 
tor in the limit e — * 0. 

The lattice displacement Aa; can be expanded in terms 
of the four basis vectors: 



Ax" 



7 M 



(16) 



hence the displacement is determined by a unique set of 
four integers TV- 5 = Ar'/e. The constraint that a path 
connects the two points of interest can be incorporated 
as four Kronecker deltas, which we express in a Fourier 
representation: 



]JS{N j ,Ax j /e) 



d 4 e 



w (2tt) 



The lattice propagator is given by 



K e (Ax) 



N=0 



.(20 



i e -iT. j 8i^ i h[A{e)] N , (18) 



where 



A{6) 



IT.*- 



JO, 



(19) 



The sum over TV of [^4(0)]^ produces every possible se- 
quence of projection operators, each with the appropriate 
exponential factor encoding the number of steps in each 
direction. When the integrals over 9j are carried out, 
only those step sequences that produce the displacement 
Ax M will survive. In particular, only the value of N equal 
to the total number of steps contributes. 

As the step size e goes to zero, the number of steps N 
for a fixed time interval goes to infinity as At/e. Con- 
vergence thus requires that the norm ||A(0)|| (i.e. the 
maximum norm of A(9) acting on a unit spinor) be less 
than or equal to unity. It is shown in the appendix that 
||j4(0)|| < 1 except when at least three of the dj coincide. 
Thus, as N becomes larger, [^4(0)]-™ converges to zero 
pointwise except at these degenerate 9j values. (Had we 
assumed an arbitrary step-speed a, convergence would 
have required a > 3.) 

When three 0j's coincide A(0) takes the form A(9) = 
\{e 101 + e l62 )Pi + e i02 P 1 , where A is the projector or- 
thogonal to Pi. It is easily seen that the action of A(8) 
then decreases the norm of any spinor except an eigen- 
spinor of Pi. [In the exceptional case where all four 9 t 
coincide A{9) is just e 10 times the identity. This leads 
to no difficulty, as the integral over 9 = 9i in l(T%|) 
just produces a Kroneker delta which sets N equal to 
the total number of steps.] Fermion doubling would oc- 
cur if Eq. i|l(J|) had an extra solution with zero frequency. 
This would be an eigenvector of A(9) with unit eigen- 
value, where 9i = —ek^n^ and 9i = 0. Inspection of 
the above degenerate form of A{9) reveals that such an 
eigenvector exists only if all 9i vanish, i.e. if fc M = 0. 

We next introduce the new variables kj := 6j/e, in 
terms of which i|18|) takes the form 



K e (Ax) = e*Y, 



N=0 



* /e d 4 k 



-Tr/e 



27T 



~A(ek)y\ (20) 



Taking the limit e — > 0, with the time interval At = Ne 
fixed, we have 



[A(ek)f = 



Eh 



N 



(21) 



= e lN ^o k i p i' 2 +0{Ne 2 ). (22) 
Moreover the limits of integration in (|20|l approach ±00, 
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hence l{20j) yields 



K e ^ (Ax) 



^ 4 E 



d 4 k 



2tt 



1 g-iTtjkjAx 3 e iNej: :i k j P j /2 



(23) 

The components Ax- 7 are defined in (|16|l relative to the 
tetrahedral basis of 4- vectors . The quantities kj can 
be viewed as components of a co- vector in the dual basis, 
and the corresponding components in an arbitrary basis 
are denoted k^, 



k^. 



(24) 



Hence we have J^. kjAx^ — k^Ax^. Substituting l|24|) 



for kj, and using (|SJ, the sum in the last exponent of 
(I23|l becomes fc M <7 M . Changing integration variables from 
kj to k^ in (|2*3f gives rise to a Jacobian |<9/c,-/<9fc M | = 
\n>j\, which can be computed using an explicit form of 
the tetrad of unit 3- vectors, yielding 



&% = 48 V3 dV. 



(25) 



The final step in taking the limit is to replace the dis- 
crete variable N by a continuous one s = Ne, in terms of 
which the sum J2n becomes J ds/e. With this replace- 
ment, and the change of variables from kj to k^, l|23|) 
becomes 



unitarity are necessarily in conflict. Indeed, as shown 
in 0, one can write a unitary discrete evolution rule 
on a body centered cubic lattice (or on an alternating 
pair of simple cubic lattices) whose continuum limit is 
the Weyl equation. (Interestingly, unitarity and locality 
were shown there to imply the Weyl equation.) Consider 
an initial state that is non-vanishing only at one lattice 
point, with normalized spin state \ip). At the next time 
step according to 10 it has support at the four corners 
of a tetrahedron, with the amplitudes ^Pi\^}. The norm 
of the state after one step is then £\ jPi\ifj) — s, i.e. 
it has decreased by a factor of two, violating unitarity. 

It is not just the norm change that violates unitarity. 
Also, the evolutions of orthogonal states do not remain 
orthogonal. Two points at one time have either one or 
no common points in their one-step future. In the for- 
mer case, the one-step evolutions of two orthogonal states 
concentrated on the two initial points are clearly not or- 
thogonal, because they overlap in just one point which 
will make the unique non-zero (since PiPj 7^ 0) contri- 
bution to the inner product. The evolutions therefore 
have "more overlap than they should" , which presumably 
counteracts the loss of norm of each individual evolution 
in such a way that the continuum limit is unitary. 
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K^ {Ax) 



48v / 3e 3 



ds 



jPk_ 



(26) 

Except for the peculiar factor in front, this is just the con- 
tinuum retarded propagator, albeit in a perhaps slightly 
unfamiliar form. The integral over fco produces a Dirac 
delta function S(s — Ax°), after which the s-integral sets s 
equal to Aa; (assuming Ax° > 0), yielding the retarded 
propagator in the more common guise of a three dimen- 
sional Fourier transform. Had we kept the subleading 
terms of order A^e 2 (= se) in l|22|) the convergence factor 
for the integration limit s — > co would have been supplied 
much as in jjj. 

It remains to account for the prefactor 48-\/3e 3 . We 
computed the propagator to go between two points on 
the hyperdiamond lattice. In the continuum, the ampli- 
tude to arrive at one point starting from another point is 
zero, since only by integrating over a finite region should 
a nonzero amplitude arise. The prefactor is none other 
than the volume per point in the lattice, with the step 
length a equal to 3e. Hence what we have actually ob- 
tained is the continuum propagator integrated over the 
volume associated with the initial lattice point. 

We have shown that a very simple discrete path in- 
tegral converges in the continuum limit to the retarded 
propagator for the Weyl equation. Although the under- 
lying lattice is not Lorentz invariant, that symmetry is 
recovered by the propagator in the continuum limit. 

Another symmetry not possessed by the discrete prop- 
agator is unitarity. This is not because discreteness and 
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APPENDIX 

In this appendix we prove that the norm of the matrix 
A defined in Eq. i|19|) is l ess than unity unless at least 
three of the 9i coincide, in which case the norm is unity. 
The proof is due to Eli Hawkins. 

Let \v) be any unit spinor. The squared norm of A\v) 
is ||A|z/)|| 2 = (u\ A'A\u), whose maximum is the larger 
eigenvalue of A^A. This value defines the squared norm 
\\M 2 . 

Using the definition of the spin projection operators 
(|u]l and the inner products of the unit vectors hi (1 if 



i = j and — ^ if i ^ j) we find 



tr(AU) = 1 + 12 ' 

(y) 



(A.l) 



where the sum is over the 6 choices of {i, j} C {1, 2, 3, 4}. 
This trace is at most 2, and therefore the smaller eigen- 
value of A* A is less than 1 unless A' A = 1. Hence 



$ := det (A* A - 1) 



(A.2) 



has the same sign as 1 — \\A\\ 
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The matrix A* A is linear in terms of e*( 9i ~ ej \ therefore 
<1> is quadratic. Because $ is invariant under all permuta- 
tions of the #'s, it can be written as a quadratic function 
of the cosines cos(#i — 6j). Because $ vanishes when the 
#'s are all equal, it is convenient to write it in terms of 
the cosines minus 1. It thus takes the form, 

$ = a^(l-cos[e». i -6» J ])+6^(l-cos[6l J -^]) 2 

+cJ2 (1 - «w[0_ - 9j\) (1 - cos[fl fc - 9{\) .(A.3) 

(ij)(ki) 

The last sum is over the 3 partitions of {1,2,3,4} into 
pairs. By setting 9\ = 02 — #3 in this expression, we see 



that a — b = 0. To determine the value of c, consider the 
case that 9i = 9 2 = and 63 = 64 = tt. Then A is the 
hermitian matrix \{fi\ + h-z) -a, which has eigenvalues 
±1/V3, hence $ = 4/9. The last sum in l|A.3|l is 8, so 
c = 1/18. 

The determinant <!> is thus given by 
$ = S £ (l-cos[^-^])(l-cos[^ 

(ij){kl) 

(A.4) 

This satisfies $ > 0, therefore ||yl|| < 1. Each term is 
non-negative, therefore $ = only if every term vanishes. 
This occurs only if at least three of the 9's are equal. 
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